Low energy electron-phonon effective action from symmetry analysis 
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Based on a detailed symmetry analysis, we state the general rules to build up the effective low 
energy field theory describing a system of electrons weakly interacting with the lattice degrees of 
freedom. The basic elements in our construction are what we call the "memory tensors" , that keep 
track of the microscopic discrete symmetries into the coarse-grained action. The present approach 
can be applied to lattice systems in arbitrary dimensions and in a systematic way to any desired 
order in derivatives. We apply the method to the honeycomb lattice and re-obtain the by now well- 
known effective action of Dirac fermions coupled to fictitious gauge fields. As a second example, we 
derive the effective action for electrons in the kagome lattice, where our approach allows to obtain 
in a simple way the low energy electron-phonon coupling terms. 



PACS numbers: ll.lO.-z, 11.30.-j, 71.10.Fd, 63.20.kd 
I. INTRODUCTION 

The synthesization of graphene has triggered an enor- 
mous interest in the analysis of its low energy properties 
using techniques from relativistic field theory due to the 
fact that electrons in this material show a Dirac like be- 
havior for stoichiometric fillings [1 . The origin of this lies 
in the band structure of fermions in the honeycomb lat- 
tice: the Fermi surface degenerates at a number of points, 
known as "Dirac points", close to which the dispersion 
relation takes a conical shape. 

The conical shape of the Fermi surface is not exclusive 
of graphene, but it shows up in a number of unrelated 
cases: fermions on the kagome lattice close to 1/3 filling 
(or 2/3 depending on the sign of the hopping amplitude) 
[H [3] and fermions in the 7r-flux state on the square lat- 
tice [4] are just two examples of a long list. One also 
encounters a degenerate Fermi surface in the mean field 
description of strongly correlated systems for particular 
solutions of the self consistent mean field equations, or in 
the Schwinger fermions mean field theory of Heisenberg 
models [5]. Last but not least, graphyne arises as another 
candidate to show up this kind physics [6, 7 . In ref.fS^ 
the appearance of Dirac cones on the Fermi surface was 
investigated in a general way, while in [9] their existence 
was exhaustively studied for all possible non-magnetic 
3-dimensional crystals. We note in passing that Dirac 
cones were observed at interfaces of topological insula- 
tors Hoi [11]. 

In the study of graphene, different approaches have 
been employed to construct the low energy effective dy- 
namics for fermions in the presence of phonons. On the 



one hand, one may start from the microscopic model and 
follow, for the case in which the lattice is deformed, the 
same steps that lead to the Dirac like Hamiltonian in 
the undeformed lattice [12] [13] . We call this the tight- 
binding approach, see [T for a detailed review and an 
extensive list of references. On the other hand, one 
may promote the configuration space in which the Dirac 
fermion propagates into a curved manifold, and then use 
the standard coupling of relativistic fermions to curved 
backgrounds [T4[ FlB j . We call this point of view the covari- 
ant approach, see the review [16^ and references therein. 
This last approach allows the description of topological 
defects jT^ as well as smooth deformed manifolds |18j[T9]. 
Despite the strong bases of both approaches, there are 
some apparent discrepancies between them, a step to- 
wards solving which can be found in \W . See also refs. 
[20 for related work applying the covariant approach. 

A crucial ingredient for the construction of the low en- 
ergy effective action is the symmetry group of the under- 
lying microscopic theory. A symmetry based approach 
has been applied to graphene has been proposed in ^T] 
and used in several other works [22H25] . The symmetries 
of a lattice system are given by the discrete translations of 
the Bravais lattice, the discrete rotations and reflections 
known as point symmetries, and possible internal sym- 
metries such as spin or color. In the process of taking the 
continuum limit to describe the physics at sufficiently low 
energies, one may naively expect that any fingerprint of 
the underlying discrete space would be washed out, leav- 
ing a quantum field theory with no detailed information 
of the particular system at hand. In other words, one 
may expect the symmetries in the continuum limit to be 
given by continuous rotations and translations, together 
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with the internal symmetries. 

In the present paper we revisit in detail how symme- 
tries are used to constrain the low energy dynamics of 
a generic lattice system, under the assumption that one 
needs to retain specific information about the underlying 
discreteness of the lattice symmetries. The main idea is 
that in the infrared, the discrete translations of the mi- 
croscopic Br avals lattice become continuous translations, 
while the point group retains its discrete nature. In other 
words, the low energy dynamics will be described by a 
field theory with continuous fields transforming in repre- 
sentations of the discrete point group. Since any term in 
the Lagrangian should be invariant under the full symme- 
try group, the system recalls its discrete origin through 
the invariant tensors of the point group that are needed 
in order to contract the fields indices. We call such in- 
variant tensors the "memory tensors" of the lattice. In 
this way, a low energy description of a system microscop- 
ically defined by certain degrees of freedom in the square 
lattice (whose point group is D4) or in the hexagonal lat- 
tice (whose point group is Dg) will differ in the content of 
memory tensors and this would lead to different invariant 
local terms in the effective Lagrangian. 

The paper is organized as follows: In Section [IT] we 
give a general description of the methodology, in Section 
|III| we show the method at work in the case of the De 
point group, and we apply it to analyze both graphene 
and kagome systems. We study as a first example the 
graphene case, and show how to derive the low energy 
effective field theory reproducing the tight-binding ex- 
pressions [12 , 13 . These results have also been derived in 
[2T] following a related symmetry based analysis. The ad- 
vantage of the present approach is that it can be applied 
to study more general cases in a straightforward manner. 
In a second step, we apply the method to fermions in the 
kagome lattice, which describes a tight-binding model at 
filling 1/3 (or 2/3, depending on the sign of the hopping 
amplitude) [2j|3]. In this example, the method shows its 
power since the derivation of the effective action is quite 
straightforward. Section [IV| contains the conclusions and 
outlook and in the Appendix we present an alternative 
way to derive the transformation properties of the low 
energy fields. 



II. GENERAL CONSTRUCTION 



A. Point group field theory 

A d-dimensional lattice system is defined by a unit cell 
containing A atoms or molecules at positions {ra}ae[i---A] 
within the cell, which is repeated periodically on a Bra- 
vais lattice, consisting in linear combinations of a set 
of d basis vectors - with integer coefficients, 

R = Yli=i ^i^i with rii G Z. In other words, such a 
system has a configuration space consisting on sites 

^Ra = ^^^a. (1) 

where the components of the vectors fa in the -(i] 
basis are smaller than one and positive. For the kind of 
lattices we are interested ir|^ a space symmetry transfor- 
mation acts on the lattice sites as 



where t is an arbitrary Bravais lattice vector, and P be- 
longs to the point group Gp which is a finite subgroup 
of the (i-dimensional orthogonal group 0{d), consisting 
of discrete rotations and refiections. The group of space 
symmetries can be decomposed as Gt >^Gp^ where Gt is 
the group of discrete translations of the Bravais lattice. 

1. Low energy symmetries and field theory 

In the continuum limit, the Bravais sites R get replaced 
by continuum Cartesian coordinates x in a d-dimensional 
Euclidean space R^. In other words, the group of trans- 
lations Gt is replaced by continuous translations R^. 
This implies that the low energy degrees of freedom are 
well described by fields ^{x) which are functions on R'^. 
Naively one would expect that the point group would be 
replaced by the continuous Euclidean point group 0{d), 
resulting in a space symmetry group R^ x 0{d). If this 
expectation were realized, the fields would transform un- 
der 0{d) and any memory of the ultraviolet discreteness 
would be lost in the infrarecQ 

The point of view taken in this paper is that the last 
step in the above reasoning must be omitted, the discrete 
group Gp must be maintained in the low energy limit as 
the group of point transformations. In other words, our 
guiding principle is that the low energy space symme- 
try group of the theory is R^ x Gp. The low energy 
degrees of freedom will therefore be described by fields 



In this section we present the procedure to con- 
struct the low energy effective action of an arbitrary d- 
dimensional lattice system. We give the guidelines to 
build invariants out of the fields representing the low en- 
ergy degrees of freedom and the memory tensors that 
keep track of the discrete point group symmetry in the 
infrared description, and discuss the basic elements defin- 
ing the field theory for electrons and phonons on a lattice. 



^ We concentrate on "symmorphic" space groups, which corre- 
spond to the case where the subset of discrete rotations and 
reflections is a subgroup of the space group. In the non- 
symmorphic case, an additional translation up that depends on 
the transformation P should be added to the right hand side of 
eq. ([2} in order to map the lattice into itself 26 . 

^ In these considerations, the group 0(d) should be replaced by 
SO(d) if the low energy dynamics is not parity invariant. 
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^£{x) with i e [1 • • • L], transforming under a point group 
transformation P G Gp according to an L dimensional 
representation of Gp 



(3) 



The low energy dynamics of the system has to be in- 
variant under xi Gp. Invariance under translations is 
guaranteed if any term in the Lagrangian is written as 
an integral of local products of the fields ^i{x) and their 
derivatives. In order to construct invariants under Gp we 
need to contract the i indices with suitable Gp invariant 
tensors, which we shall call "memory tensors" . A generic 
Gp-invariant tensor of order p satisfies 

Te,e,...i, = ^ [D{P)]i,iAD{P)]i,i, • • • [D{P)U_^i,Te^e,...e^ . 



(4) 

Thus a generic x Gp invariant term in the Lagrangian 
density, of order p in the fields and containing no deriva- 
tives, reads 



Cqo 



E 



Ti,i,...i^^i,{x)^e,{x) • • • ^i^{x) . (5) 



Notice that derivatives di transform under the vector rep- 
resentation of Gp, that in what follows we will denote as 
V. The vector representation is defined by 



(6) 



where [P]ij is the standard d-dimensional matrix repre- 
sentation of P. In consequence, terms containing deriva- 
tives can be built in a similar manner, with the help of 
invariant tensors with as many additional vector indices 
as derivatives are present. Explicitly, a term of degree q 
in derivatives and order p in the fields reads 



t^di^...^t,{x)---d...i^^i>^{x) 



(7) 

where the tensor Ti^...i^£^...£^ satisfies Q with the i in- 
dices transforming in the vector representation!^ 

2. Obtention of the invariant tensors 

The L-dimensional representation D{P) carried by the 
fields ^£{x) is in general reducible. A reducible represen- 
tation D{P) can be decomposed as a sum of irreducible 
representations of Gp as 



(8) 



^ If internal symmetries (such as spin or color) are present, the 
fields have additional indices transforming in a representation of 
the internal group, that must be contracted with the correspond- 
ing invariant tensors. 



here J labels the irreducible representations of Gp and 
the integer aj accounts for the multiplicity of the irre- 
ducible representation J in the decomposition. Among 
the irreducible representations of Gp that will be of in- 
terest in what follows we have the singlet representation 
E in which any point group transformation acts as the 
identity, and the vector representation V that we defined 
above (|6|. 

A given term in the action contains a product of fields 
and their derivatives, that transform in the corresponding 
tensor product representation of Gp. The existence of 
an invariant tensor to contract all the indices in such a 
product can be traced back to the presence of the singlet 
representation in the corresponding decomposition ([8|. 

To obtain the decomposition of a given representa- 
tion D{P), one needs to compute the character vector 
= (tr(L>(Pi)), tr(L>(P2)), • • • where Pi, P2 are repre- 
sentatives of the conjugacy classes Ci, C2, the dots stand- 
ing for the remaining classes of the point group Gp. The 
multiplicities aj are obtained by linearly decomposing 
in terms of the characters of the irreducible represen- 
tations 



n 



E 



(9) 



Given a representation D{P), one finds its projection 
onto the irreducible representation J by use of the pro- 
jector |37] 



dim J 

TgT 



^ tr{D'{P)fD{P). 



(10) 



PeGp 



In order to construct an invariant tensor as those entering 
into eqs. ([5| and ([7|), one needs to project the tensor 
product representation into its singlet component, using 



1 



\Gp\ 



(11) 



PeGp 



The action of this projector is easy to understand: acting 
on an object that transforms under D{P)^ it gives the 
sum of all its images under point group transformations. 
The resulting object will necessarily be invariant under 
Gp. Explicit examples will be given below. 



B. Phononic and electronic fields 

The natural question is then which representations of 
G p have to be considered when constructing the effective 
field theory of a given lattice system, or in other words, 
what is the correct field content ^^(x). In this section 
we obtain the field content that takes into account the 
low energy degrees of freedom related to phonons and 
electrons moving in the lattice. 
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1. Phononic field 



2. Electronic field 



To analyze the vibrational degrees of freedom of the 
lattice we denote by via{R) the displacement of the 
atom a on the Bravais site R in the I-th direction of 
3-dimensional space / G [1---3]. Notice that for the 
case of a lower dimensional d < 3 lattice, the vja include 
the transverse displacements as well as the displacements 
along the directions of the lattice. They give a total of 
3A degrees of freedom which enter quadratically in the 
Lagrangian at the harmonic level. The resulting 3A nor- 
mal modes can be classified, according to the standard 
theory of phonons in solids, as 3A — 3 optical modes plus 
3 acoustic modes [27 . At sufficiently low energy only the 
acoustic modes are relevant, and they have linear disper- 
sion relation. In consequence, at low energy we need to 
consider only three fields vi{x). 

The fields vi{x) can be decomposed as vj = 
(lii, • • • ^Ud^hi^ - • • , hs-d). The Ui{x) fields take into ac- 
count the deformations of the Bravais lattice and trans- 
form under the point group Gp in the vector representa- 
tion V defined in (|6|. On the other hand, the 3 — d fields 
hr{x) account for the transverse (or "fiexural") displace- 
ments, and they are scalars i.e. transform in the singlet 
representation E. 

To isolate the elastic degrees of freedom from the global 
translations and rotations of the lattice on the Euclidean 
3-dimensional space in which it is embedded, the effective 

action must depend only on the combinations dihr{x) 

and Uij{x) = l/2{diUj{x) djUi{x)) [28^. The field P acts on Bloch states as 



In order to derive the transformation properties for the 
electronic degrees of freedom, we need to work out the 
implications of eq. ([2|. To do that, we note that under 
a point group transformation P, the position fa of the 
a-th atom inside the unit cell transforms as 



(17) 



Then under a point group transformation the atom a in 
the Bravais site R is mapped into the atom ap in the 
Bravais site P • R -\- . The Bravais lattice vector 
arises whenever the action of P takes fa out of the unit 
cell (see Fig. [T|). 

Calling Ca{R) the operator that annihilates an electron 
in the Wannier state localized around the site x^^, a 
point group transformation P acting on the lattice maps 
such state into the corresponding Wannier state centered 
around the transformed point P • ^p^, namely]^ 



c'aiR) 



Bloch states are annihilated by the operators 



e-*^Ca(^). 



(19) 



(20) 



Bravais 



Using (19), one finds that a point group transformation 



dihr(x) transform in the vector representation 
diK{x) = Y,Piidihr{P~^-^) 



and we denote this as 



dih e V . 



(12) 



(13) 



c'aik) = E 

i?G Bravais 

In other words we obtain 



CaAR)- (21) 



c;(fc) = e*(^-'=)-*« c„,(P-fc). 



(22) 



Regarding Uij{x)^ it is a symmetric tensor that trans- 
forms in the product representation V (g) V as 



As will become clear below, the phase factor on the right 
hand side characterizes the transformation properties of 



in other words 



Uij e S{\ (g) V) , 



(14) 



(15) 



where S stands for "symmetric part". Notice that any 
symmetric tensor can be decomposed into its trace and 
its symmetric traceless parts, namely 



= ^iX/e/c(^) + ^ (^ik^jl - 

k kl ^ ' 



Ukl{x) . 

(16) 

Since the first term is invariant under point group trans- 
formations, the representation 5(V V) contains a sin- 
glet in its decomposition, i.e. <S(V (g) V) = E • • • . This 
singlet describes the homogeneous dilation of the lattice. 



^ To prove this, let us write the Wannier wavefunction centered 
around the atom a on the site R as = '^a(f — R — ra)^ an 

electron in this state is annihilated by the operator Ca{R)- Ap- 
plying a point group transformation P to the space coordinates 
r we get 

= lPa{P~^ -f- R-ra) = 

= i^aiP-^ ■ (f- P ■ R - P ■ fa)) = 



■ i'apif- P R- Tap 



-tj)- 



(18) 



where in going from the third to the fourth line, we used ( |17| and 
the fact that the function ipaif) satisfies ipap(f) = ipa(P~ • f)- 
This implies that an electron in the transformed wavefunction 
will be annihilated by Cap{P • R-\-t^). 
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the low energy electronic degrees of freedom under point 
group transformations. 

In a non-interacting system, the Hamiltonian is diago- 
nalized by linear combinations of Bloch states with fixed 
k. Since the index a takes A values (corresponding to A 
atoms or molecules in the unit cell), there would generi- 
cally be A different energy bands {Sn{k)}ne[i---A]' Calling 
Cn{k) the operator that creates an energy eigenstate with 
crystalline momentum k in the n-th band, we can write 



Ca(fc) = Yl 



tta 



(23) 



n=l 



Bloch theorem implies that states with crystalline mo- 
mentum k accommodate in representations of the group 
Gj^ consisting in those point group transformations that 

leave k invariant, the so called "little group" of ^0 This 
in turn implies that, if there exist values of k such that 
the transformation rule (22) contains a non-trivial repre- 



sentation of the little group G^, then the n indices mix 
under little group transformations, and the energy bands 
will be degenerate at those points. We will assume that 
such degeneracy occurs at isolated points {km}me[i---M]- 

The ground state of the non-interacting electron sys- 
tem is obtained by filling all Bloch states with nega- 
tive energy en{k) < while those with positive energy 
en{k) > are empty. We will assume that interactions 
do not alter this picture. 

Now we are ready to take the low energy limit. Given 
an energy cutoff £ cutoffs high energy bands, those which 
lie completely above it, are not excited at all at low ener- 
gies. On the other hand, low energy bands can be excited 
only for momenta satisfying |£n(^)| < ^cutoff- 

We will assume in what follows that low energy bands 
have their minima or maxima at degenerate points 
{^m}mG[i -M] ^ud wc will tuuc the chemical potential 
to ensure that those points lie at zero energy. In conse- 
quence low energy states have momentum close enough 
to some of the krn- To construct the low energy effective 
action we replace the annihilation operators by 



^nra 



(24) 



2i k-V - 



^ The proof is as follows: Bloch's theorem states that the eigen- 
functions of an electron in a Bravais lattice can be written as 
'0^g(r) = exp(zA: • r)w^g(r) where the function w^g(r) is periodic 
in the lattice and satisfies 

/_ 
.2m V 

This equation must be solved in the unit cell with periodic 
boundary conditions. Although the potential V(f) is invariant 
under the action of the full point group Gp, the equation is only 
invariant under the restricted set of point transformations that 
satisfy P • k = k, i.e. under the little group of the momen- 
tum k. The eigenfunctions '^^^(f) therefore accommodate in a 
representation of G^. 



With these operators we can recover, in our original basis 
([23) 



an ^nm 



(25) 



Notice that in (24) we have also expanded the high 
energy bands annihilation operators around the minima 
of the low energy bands, and that we have included the 
resulting operators in the sum of (25). Since this cor- 
responds to an assumption on the high energy modes, 
it would be harmless as long as the corresponding ad- 
ditional degrees of freedom are gapped in the effective 
theory. This observation will be important in Section III 
when we study the kagome lattice. 

Since under point group transformations the energy 
bands have to be invariant, then the transformation P of 
one of the minima km must necessarily result in another 
minimum kmp^ namely 



P ' km. — kr) 



(26) 



With this, and using formula (22), we can state that un- 
der point group transformations low energy annihilation 
operators transform according to 



(27) 



where, in the low energy limit, we have disregarded the 
q dependence of the exponent in the right hand side of 
(22). By Fourier transforming this low energy operators, 



we get our electronic annihilation fields 



*ma(^) 



(28) 



The coordinate x in the above equation denotes the posi- 
tion in the Bravais lattice. The transformation rule (27) 



determines the representation of the point group under 
which the electronic fields transform, its properties being 
fixed by the transformation rule of the atom positions fa 
and of the minima krn of the low energy bands. Explicitly 

*:„a(^) = / d'qe^'^-^c'^M 

= e'^'"p<"«'^,a,(P-f) (29) 

To make the above equation compatible with the general 
transformation rule ([3| of a field, we need to replace P 
in the above formula with P~^, to obtairj^ 



(30) 



^ This is consistent with the fact that transformation on the 
wavefunction coordinates r acts as a transformation P on the 
lattice. 
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with 



A. De group and its representations 



[D(P-I)] 



ma ma — ^mm p ^aa p ^ 



(31) 



where the indices m, a determine mp^ap. We denote this 
representation as R, i.e 

^ma^R. (32) 

Notice that this representation is in principle reducible. 

To summarize, in order to construct the low en- 
ergy field theory we need to identify the discrete point 
group Gp. Low energy fields representing phonons con- 
sist on fiexural dihr{x) and in-plane Uijix) degrees of 
freedom, transforming in the vector V and symmetric 
product 5 (V V) representations of Gp respectively, 
while those representing electrons consist on ^rna(^) 
transforming in the representation R defined in (30). 



In consequence our field content is given as ^i{x) 
{dihr{x)^Uij{x)^^ rna{^))i ciud invariant terms will be 
built out of memory tensors with indices in the corre- 
sponding representations. The list of such tensors can be 
exhausted by group theory methods, thus providing the 
most general effective action compatible with the point 
symmetries. Up to this point we have only made use of 
the point symmetry group of the system. At low energies, 
the translation subgroup enters only through its contin- 
uum description, which implies the conservation of q. 
Nevertheless, the underlying discrete translation symme- 
try imposes the conservation of crystalline momentum, 
thus a term in the Lagrangian originating transitions be- 
tween states at different minima km in the Brillouin zone 
must be forbidden. This implies that, among all the pos- 
sible invariant tensors, we must keep only those which 
are completely diagonal in m. 

In the next section we apply this procedure to a par- 
ticular important case. 



III. De FIELD THEORY 

In this section we apply the methodology described 
above to the case in which the point group is Dg, which 
arises, in particular, for the cases of the hexagonal and 
kagome lattices. 

We first discuss the common features of these two dif- 
ferent lattices, i.e. the decomposition of reducible rep- 
resentations into irreducible ones that allows us to con- 
struct the invariant tensors. 

In a second step we construct the singlet terms (Dg- 
invariants) that conform the low energy effective action 
as an expansion in momenta. The results are compared 
with the literature. At this points we analyze the differ- 
ent lattices separately, showing how to derive the relevant 
invariant terms in each case. 



In general, the "dihedral group" is defined as the 
symmetry group of a regular n-polygon. It is generated 
by two elements R.,Fx as 



{R,F, I = e, = e, F'^RF, = R'^} 



(33) 

here e is the neutral element of the group. 

The Dg group is a non-abelian discrete group of order 
I Del = 12. It can be generated from a 60^ rotation R and 
a refiection through the x-axis Fx. It has six conjugacy 
classes: those of the identity and R^ of dimension 1, those 
or R and R^ of dimension 2 and those of F^ and RF^ of 
dimension 3. 

For completeness we mention that De is a product 
group Dg = D3 X Z2 with D3 the symmetry group of an 
equilateral triangle generated from a 120^ rotation and a 
refiection. 



1. Representations 

The De group has six irreducible representations, that 
we denote by 



J = E,Ai,A2,A3,V,V\ 



(34) 



where E is the singlet representation, Ai, A2, A3 are 1- 
dimensional and V,V' are 2-dimensional. V is the vec- 
tor representation (|6| and we refer to Ai and V as 
"pseudoscalar" and "pseudovector" representations re- 
spectively. 

The pseudoscalar representation Ai, which will appear 
below, corresponds to representing the rotations trivially 
and the refiections by a sign. 

The (faithful) vector representation V is written natu- 
rally in terms of 2 x 2 orthogonal matrices acting on the 
plane 



/ 1 _Vs 



D^{Fx 



2 2 
1 
-1 



(35) 



(36) 



The (non-faithful) pseudovector representation V can 
be obtained from the V matrices as 



D^\r) = {D^{R))^ =12 2 



^/3 _ 1 
2 2 



Due to non-faithfulness of V one has 

= {D'^'{R)fD^'{F,) = 
= D^'iF,) , 



(37) 
(38) 



(39) 
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where we have used (37) to go from the second to the 
third hne in (|39|. Hence, reflections through x and y- 
axis act on a V^^pseudovector in the same way. 

In what fohows, to distinguish vectors in the 2- 
dimensional representation V from pseudovectors in the 
2-dimensional representation V', we use for the first 
Latin indices, while for the second we use Greek indices. 

The decomposition of a given De representation D{P) 
into irreducible representations amounts to compute the 
character 



X 



D 



(tr(e), tr(i?3), tr(i?), tr(i?2), tr(F,), tr(i?F,)) , (40) 



where tr(P) = tr(I)(P)) with P G De, and to write it as 
a linear combination of the characters of the irreducible 
representations. The character table for Dg can be found 
in Ml 





Ce 


Cr3 


2Cr 


2Ci^2 


3Cf. 




E 


1 


1 


1 


1 


1 


1 


Ai 


1 


1 


1 


1 


-1 


-1 


A2 


1 


-1 


-1 


1 


-1 


1 


As 


1 


-1 


-1 


1 


1 


-1 




2 


2 


-1 


-1 








V 


2 


-2 


1 


-1 









The horizontal lines of this table contain six charac- 
ter vectors fo^" each irreducible representation. 
These vectors are orthogonal with the metric gij = 
diag(l, 1, 2, 2, 3, 3) whose diagonal entries are given by 
the dimension of each of the classes (quoted in the top 
line of the table). 

Decomposing (40) in the basis given by {x*^} as in ([9| 



one obtains the coefficient aj giving the number of times 
the irreducible representation J appears in the represen- 
tation D{P). 

For the construction of invariants it is useful to quote 
the decomposition of the product of any two irreducible 
representations. One has 



E (g) Irrep = Irrep , 

Ai(g)Ai=E {similar for A2, As), Ai (g) V 

Ai (g) A2 = A3 {and cyclics), Ai (g V = 

A2(g)V = V^ A2(g)V 

A.sg)V = V, A.g) V 



Vg)V =E©Al©V^ 
V©V = A2©A3©V, 
V'^V = E©Al©V^ 



V, 

= v, 

= V, 



(41) 



These rules allow us to decompose and arbitrary product 
of representations. 



2. Memory tensors 

As mentioned above, the memory tensors have their 
origin in the decomposition of products of representa- 
tions into irreducible representations. As an example, the 
first line in (41) indicates that a product of two vector 



representations decomposes as a sum of three irreducible 
representations: the first two are 1-dimensional and the 
third one is 2-dimensional. 

The presence of E in the decomposition of V V in- 
dicates the existence of a tensor with which we can con- 
struct a scalar out of two vectors. To obtain such tensor, 
we apply the projector Pe to an object Oij transforming 
in the product representation, to obtain 



12 
1 



(oii + 022)^^^ 



(42) 



Stripping the scalar factor, we extract our invariant ten- 
sor in the V © V representation. In a similar way, one 
can extract the pseudoscalar part Ai as 



12 ^ 

p 

1 



s{P)PuPjjO,^ 



(Oii - 022)ei 



(43) 



where s{P) is the sign representing the transformation 
P in the pseudoscalar representation Ai. The resulting 
tensor Cij gives a pseudoscalar when contracted with two 
vectors. This tensor can be also understand as the singlet 
in the decomposition 



Ai (g) V (g) V = E © Ai © V' 



(44) 



Finally, the presence of V in the V©V decomposition 
calls for the existence of a Dg-invariant three index tensor 
Saij such that when contracted with two vectors gives a 
pseudovector. Its existence can also be inferred from the 
singlet in the decomposition 



V © V © V = E © Ai © 3V , 



(45) 



Such tensor can be obtained by using the corresponding 
projector, and it reads 



S 



lij 



{(yi)ij . 



(46) 



where (11,(73 are the standard Pauli matrices. According 
to (|4| the invariance of the tensor Saij under De implies 
that it satisfies 

Sa^j = Y}^^' \D^iP)h [D^{P)],j Sc^: ■ (47) 



This tensor turns out to be a key ingredient in the con- 
struction of the effective action. 

For completeness we observe that from the last line 
in (41) we can infer the existence of a tensor 5ocj3 that 
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provides a scalar out of two pseudovectors, a tensor CafB 
providing a pseudoscalar, and a three index tensor Tq,/3^ 
providing a pseudovector. This last tensor can be under- 
stood as the singlet in the decomposition 

V V V = E © Ai © 3V\ (48) 

Its explicit components are given by 

Tia(3 = {o-s)a(3, ^2a/3 = —{cri)a(3 • (49) 



The Tai3j tensor (49) is totally symmetric in its indices 
and traceless on any two indices. It satisfies 

(50) 



This tensor allows to obtain the singlet in (48) out of the 
product of three V pseudovectors. 



B. From memory tensors to the effective action 

Let us first analyze the symmetry properties of the low 
energy degrees of freedom we want to describe, that is 
find the field content describing electrons and phonons on 
our De field theory. First we discuss the phonon modes, 
which have the same symmetry properties for the two 
examples we want to analyze in detail, i.e. the hexagonal 
and the kagome lattices. Then we turn into the electron 
modes, which have different transformation properties in 
each system. 



1. Phonons 

As mentioned in Section [TTj the elastic degrees of free- 
dom split into in-plane phonons Uij{x) transforming in 
the symmetric tensor representation <S(V©V) and out of 
plane fiexural phonons dihr{x) transforming in the vector 
representation V. In the present case, i = x^y and the 
index r can be omitted since there is a single transverse 
direction. 

The vector representation is generated by the matrices 



(36) and it is irreducible. Then for the fiexural phonons 



we can simply write 



(51) 



On the other hand, the tensor representation V © V 
decomposes as in the first line of (41) in terms of irre- 



ducible representations. It consists of tensors with two 
vector indices, which allows us to use the standard de- 
composition of a tensor into its trace, its antisymmetric 
and its symmetric parts, in order to get an interpretation 
of the irreducible representations entering into (41). In- 



deed, as advanced in Section [II| the singlet representation 
E corresponds to the trace part. Being two dimensional, 
the V representation corresponds to the two dimensional 



space of symmetric traceless 2x2 tensors. Finally, the 
Ai representation must then correspond to the one di- 
mensional space of antisymmetric 2x2 matrices. With 
this, we can deduce 



g5(V©V) =E© V. 



(52) 



We can arrange the two independent components of the 
symmetric traceless part of Uij{x) in a V pseudovector, 
by making use of our invariant tensor Saij as 

^cx^ i-^) ~ ^ ^ ^aij '^ij (^) — 

= {dxUx - dyUy, dj^Uy + dyUj,) . (53) 

We would like to stress that A^\x) is not a vector in 
the usual sense since it transforms in the V representa- 
tion. This object has been named pseudovector since 
under a y-axis refiection (x^y) {—x^y) transforms 
as (Ai^^(f),4''^(f)) ^ {A^a^\x),-Ai''\x)) contrary to 
naive expectations. This last transformation is noth- 
ing but the expected one as soon as one recognizes that 
{x) transforms as a V representation of Dg 

/3 



where D^\P) are as in (l37|), ([38). 



2. Electrons 

In order to obtain the representation in which the elec- 
tronic fields transform, we follow the lines sketched in 
Section |ll| The two examples we study, namely that of 
graphene and electrons on the kagome lattice, correspond 
to a triangular lattice with symmetry group Dq. The 
first Brillouin zone consists of a hexagon whose vertices 
are identified under 120^ rotations. In other words, we 
have two independent vertices km with m = +, — , which 
have a non- trivial little group. As explained in Section 
\n\ this implies that at these points k± , the energy bands 
degenerate, and our technique applies as long as the fill- 
ing selects the corresponding energy as the ground state. 
This is precisely the case for graphene at stoichiometric 
fillings and in the kagome case we adjust the chemical 
potential to have filling 1/3. 

The triangular Bravais lattice is generated as ^ = 
nitti + n2a2 with ni, n2 G Z and 



ai = (l,0). 



^2 



(cos-,sm-). 



(55) 



where for convenience we set the lattice spacing to one 
by an appropriate choice of units. The two inequivalent 
Fermi points are located at 



^4 



47r 



(56) 



a. The graphene case 

Being made out of carbon atoms arranged in a hexago- 
nal structure, graphene can be seen as a triangular lattice 
with a fundamental cell containing two atoms at positions 
ra with a = a, b. The low energy dynamics is hence de- 
scribed by two energy bands. The positions of the a and 
b atoms inside the unit cell are 



1 TT . TT 

^(cos-,sm-), 



2 TT . TT 



To analyze the low energy limit we consider the neigh- 
borhood of the Dirac points k± and write k = k± -\- q 
with 1^ small enough. We expect that only the degrees 
of freedom in such regions participate in the dynamics 
and therefore we get four low energy independent anni- 
hilation fields 



^ma (^) 



(58) 



These fields transform according to (30), (31) and to 



make the transformation rule explicit, we need to com- 
pute t^, ap and mp defined in (17), (26). In order to 
do that, let us first notice that the action of the gen- 
erator R on Br avals sites is homogeneous, while for the 
unit cell basis f^^fb one gets a non- homogeneous piece 
(see Fig. [l]) 



R-ra=fb 



- a2 — tti , 



from which we extract 



and 



bi? = a , -R 



—^R ~^R 

= — ai , = a2 — ai . 



Second, for the generator we get 
implying 



and 



aF, = b , 
bF, = a , 



In consequence from (27) and (29) we obtain 
/ uj\ 



cj* 



cj* 
\uj 



(59) 

(60) 
(61) 

(62) 

(63) 
(64) 

(65) 



J 




FIG. 1. In the honeycomb lattice, a rotation P maps the 
atoms a and b corresponding to the cell at position into 
the atoms b and a corresponding to two different cells, none 
of which sits at the rotated position P R. In order to express 
the positions of the rotated atoms in terms of the atoms cor- 
responding to the rotated cell, we need to use the vectors ta 
and tb- 



/ cj* 
uj 

V 



UJ 

cj* y 



(66) 



with 6^ = 6*3". 

The character for the De representation defined by 



(65)-(66) is 



X 



^ = (4,0,0,-2,0,0), 



(67) 



This representation is reducible and decomposes in Dg 
irreducible representations as [22, 



G R = V e V . 



(68) 



Quadratic terms: 

Having obtained the representations of De contained 
on the basic low energy fields Uij{x), dih{x) and "^rnaix)^ 
our aim now is to build an invariant effective action. 

As a first step, one might wonder whether quadratic 
fermionic terms, not including derivatives, are allowed in 
the low energy effective action. This amounts to analyze 
the De singlets in the ^J^^^^no decomposition. From (41 ) 



one immediately finds that the 16-dimensional fermion 
bilinear representation decomposes as 

^L^ma G R*0R = 2E©2Ai©2A2©2A3©2V©2V . 

(69) 



The coefficient of the singlet E in equation (69) implies 



that two possible memory tensors M, M exist. Upon 
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computation one finds 



M = I4x4 

/ l\ 

1 

1 
\1 J 



M 



(70) 



The appearance of the M tensor is not surprising, since 
the fermion representation defined by (65)-(66) is easily 
seen to be unitary. By looking at the structure of the 
M tensor, one realizes that it leads to fermion bilinears 
M"^ that mix both the Fermi points ± and the sublat- 
tices a, b. As discussed at the end of Section [n[ conser- 
vation of crystalline momentum does not allow for this 
last kind of terms in the Lagrangian, thus we must not 
consider M in the low energy limit. We conclude that, at 
the quadratic order and with no derivatives, symmetry 
allows only one term in the Lagrangian constructed from 
the trivial tensor which reads 



(71) 



(here and in what follows the ma indices are not made 
explicit whenever they are summed over). The arbitrary 
constant fi can be identified with a chemical potential. 

The tensor M also allows to write a kinetic term of the 
form 



(72) 



Next, one can explore the possible derivative terms. To 
do that, we have to decompose the product "^l^^dj'^rna 
into irreducible representations. This gives 

^ma ^V(g)R*(g)R = 2E©2Ai ©2 A2 ©2 As ©6 V©6V: 

(73) 

The presence of two singlets in the above decomposition, 
whose origin can be traced back to the presence of two V 
terms in (69), can be understood as the existence of two 



independent sets of invariant "Dirac gamma" matrices 7^ 
and 7i satisfying 



(74) 



In other words, the bilinears ^'''7^^ and ^'''7^^ con- 
structed out from those 7-matrices, transform as vectors. 
As for the tensors derived above ( 70 ) , due to momentum 



conservation only one of the sets satisfying (74) is per- 
mitted in the low energy effective action. Indeed, upon 
analyzing the structure of the gamma matrices 7, 7 one 
observes that the fermion bilinear that results from the 7 
set preserves the Fermi points while the one that results 
from the 7 set mixes the two Fermi points. Again, by 



^ In order to have gap opening terms, it is necessary to relax some 
of the assumed symmetries [33H36] 



crystalline momentum conservation we must discard the 
7i set in the low energy effective action. The relevant set 
of Hermitian 7-matrices satisfying ( 74 ) reads 

\ 



-i27r/3 



7l 



pi27r/3 



-itt/S 



t/3 



72 



-i7r/6 



pin/6 



-i57r/6 



V 



(75) 



(76) 



The reason for calling these tensors 7-matrices follows 
from their anti-commutator algebra 



(77) 



and from the form of the singlets we can construct out 
of them. 

From all these ingredients it now follows that the De- 
invariant term that can be constructed out of fermion 
bilinears and one derivative is 



(78) 



where v is an arbitrary constant. As shown below, this 
expression coincides with the derivative term originally 
obtained in the literature from a tight-binding approach. 
The virtue of our approach is that it makes manifest the 
transformation properties of the fermionic excitations un- 
der the point group. 

Accidental rotation symmetry: 



From the fermionic terms we just constructed ( 71 ),( 72 ) 



and (78), we conclude that, at the quadratic leve. 



m 



fermions, the Lagrangian for graphene reads 

£2 = ^^■^St^ -^^^Yl ^^7^^^^ + A^^^^ i^^) 

This Lagrangian can be shown to have an accidental 
continuous symmetry, already discussed in where 
it was called "full intravalley rotational symmetry". In 
our formalism, this symmetry arises naturally from the 
fermionic rotation generator A, constructed from the 7 



A = ^ [71,72] = ^ 



/l 



V 



-1 



1/ 



(80) 



From it we can define the continuous matrix D{0) as 



m 



(81) 



7- , 

e*2 y 
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which satisfies the relation 

iDiO))^^iDiO)=Ri,iO)j, 

with 



R{0) 



COS — sin 6> 
sin cos 



(82) 



(83) 



Expression (82) ensures that the quadratic Lagrangian 
£2 is invariant under the simultaneous transformations 



f R{0) - f , 



^ma(^) 



[D{0)]r 



.^^rnaiiRm-'-x). (84) 



Notice that upper and lower components of "^rnai^) do 



not mix under (81 ). 



Being an accidental symmetry of the lower derivative 
quadratic action, this symmetry does not need to be pre- 
served by interactions nor by higher derivative terms, and 
as we will see in the forthcoming sections, it is indeed 
broken by electron-phonon interactions. 



Electron-phonon couplings: 

To construct the invariant terms coupling phonons to 
electrons, let us first concentrate in the in-plane phonons 
Uij{x)^ which transforms in the 5(V(8)V) representation. 
We hence need to look for singlets in the decomposition 
of the product 



^ 5(V © V) R* ^ R = 

= 4E © 4Ai © 4A2 © 4A3 I 



(85) 

) 8V © 8V' , 



In the first line of this equation, we can use (52) and 
(|69|) to advance that we will have four singlets in the sec- 



ond line: two of them coming from the singlet in (52) 



multiplied by the two singlets in (69), and the other two 



coming from the V' in (52) multiplied by the two V in 
(69). Similarly to what happened with the quadratic in- 



variants, momentum conservation precludes half of them, 
therefore in what follows we only write down the two ten- 
sors which do not mix the Fermi points. The one coming 
from the product of the singlet in (52) with the singlet 
in (69) reads 



jT- — 



(86) 



This term has its physical origin in the changes in the 
area of the unit cell due to contractions and dilatation^ 
The singlet arising from the product of the V in 



with the V' in (69), in complete analogy with the discus- 
sion above on the derivative terms, implies the existence 



See eqs. (150) and (169) in [T] 



of a set of gamma matrices 7q, that project the fermion 
bilinear (69) onto its relevant V component. They read 



71 



t/3 



t/3 



-ZTT /3 



72 



and satisfy 



V 
/ 

V 



oi7r/3 



(87) 



t/6 



oi7r/6 



-ITT /6 



I 



{D{P))^%D{P) = D^;{P)jp. 
The resulting interaction term has the form 



(89) 



(90) 



where the pseudovector was obtained in ( 53 ) by pro- 



jecting the phononic components into the V representa- 
tion making use of the tensor Saij defined in (46). This 



interaction term reproduces the expressions previously 
obtained in the literature [1^. 

Now we can turn into the invariant terms coupling fiex- 
ural phonons to electrons. To lowest degree in fields and 
lowest orders in derivatives, we need to look for singlets 
in the decomposition of the product 



Si/l^L,^maeV©R*©R. 



(91) 



Since the irreducible representations entering into this 
product are exactly the same as in (73), we repeat the 



analysis we made to write the quadratic derivative terms 
and obtain the invariant coupling 



(92) 



To the same order in the fields, we can easily go to 
higher orders in derivatives. Since djdjh transforms in 
the same representation as Uij, eq. (85) and the subse- 



quent considerations immediately tell us that the allowed 
invariant couplings are 



and 



where 



(93) 



(94) 



(95) 
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The interaction term (94) has been recently shown to 



appear when a couphng to the spin degrees of freedom is 
present [30] . 

Notice that in the absence of a substrate, the trans- 
formation h{x) —h{x) must be a symmetry of the 
Lagrangian, which imphes that the couphng constants 

^(/i) and must be set to zero. In that case, the 
first non-trivial interaction arises to quadratic order in 
dih{x). Again since dih{x)djh{x) transform in the same 
representation as Uij {x) , we can write down the relevant 
terms without any additional analysis. They read 



and 



where 



=^>/.)^^^(V/l)^ (96) 



C,^,,,=qihh)^^^la^Ai^^\ (97) 



^i^""^ (^) = Yl ^-^J dM^)djHx) . (98) 



The procedure just presented can be pursued to any 
degree in the fields and to any order in derivatives, al- 
lowing for the construction of the most general effective 
action for electrons coupled to phonons in graphene. 

Final result and comparison with the literature: 

Then, as a general conclusion, to the lowest order in the 
fields and their derivatives, the most general Lagrangian 
consistent with the symmetries read 

jC = i^'^dt'^ ^ivY^ ^"^Ti {di - igdih) ^ + 




(99) 

The first observation regarding this Lagrangian is that 
even if the first two lines in ( 99 ) are invariant under the 
accidental continuous symmetry (84), the third line is 
not. This is due to the fact that the tensor Saij en- 
tering into the construction of the pseudovector fields 
Ao,{x) is not invariant under R{0) but only under trans- 
formations in Gp. A second observation is that the 
coupling to dih{x) in the first line can be reabsorbed 
by a local phase redefinition of the electronic fields 
^ma(^) = exp(— z^/i(x))^rna(^)- Nevertheless, such op- 
eration would give rise to a kinetic coupling —g^'^^dth. 
A final point is that, since the Ao^{x) fields are pseu- 
dovectors, a local phase redefinition of ^ma{x) cannot 
be reabsorbed by a derivative shift in any of them. In 
other words, the fields Aq,(x) are not "gauge" fields in 
the standard sense. 



We close this section showing the consistency of our 
results with the existing effective Lagrangians computed 
from a tight-binding approach, to lowest order in momen- 
tum and deformation fields. In order to make contact 
with the literature, we rewrite the electron field as 




From our kinetic Lagrangian Cy in ( 78 ) we can extract 
the quadratic part of the Hamiltonian density 

Hy = -iv Y (^+ + ^- ^idii^-^ , (101) 

i 

where = {cr 1^(12) are the usual Pauli matrices, = 
(cTi, —(J2). This coincides with the result obtained in [12]. 

As for the electron-phonon coupling, the term arising 
from (86) reads 



(102) 



and coincides exactly with the one presented in ref. 
Finally from (l90]) one obtains 



^?(.) = -i(u) E - ^-^ci'-) ^i") ,(103) 

a 

where da = (cri,cr2) and = (cri,— cr2)- This again 
coincides with eq. (3.9) of [12\. 

b. The kagome lattice case 

We now sketch the application of our technique to the 
kagome lattice, highlighting the main differences with the 
graphene case discussed above. 

The kagome lattice can be seen as a triangular lat- 
tice with three atoms per fundamental cell. As for the 
graphene case we consider the simplest non-trivial case 
of a doubly degenerate low energy band structure with 
two inequivalent Dirac points, De symmetry implies that 
they should be located at the vertices k± of the Brillouin 
zone. The third band is not degenerate and it will be 
considered a high energy band. This instance in particu- 
lar includes the description of the kagome lattice at 1/3 
filling. In consequence, we have six degrees of freedom at 
low energies, that we arrange as 



(104) 



To derive the transformation properties of the six de- 
crees of freedom, we consider the triangular lattice to be 
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generated by (55). The position of the a = a, b, c atoms 

1 



-(ai + a2), fb = 7:^2, 



-ai. 



(105) 



inside the unit ceh are 
^ _ 1 

The action of the Dq generators i?, Fx on the basis fa is 
R-fa = fc^a2-ai, 

R-fc = ft>. 



(106) 



and 



Fx-fb = fa-a2, 

Fx • fr = fr , 



(107) 



from where one reads the information needed for com- 



puting (27), as 



ai? = c, 
bi? = a, 
CR = b. 



(108) 



and 



t3^ = a2-ai, =-^1. C = 0, (109) 

for the R generator, while for the Fx generator we get 

aF, = b , = - , 

bF, = a , -F, + , 
Cf. = c , 

and 



-a2 . 



t. 



(110) 
^ =0.(111) 



In consequence we obtain 



D{R) 



cj* 0\ 
1 




and 



D{Fx) = 



uj 
1 
Vcj* 



/ 

cj* 
1 



(112) 



v 



cj* 

1/ 



(113) 



The character for the representation defined by (112) 
(Iml) is 



X 



R 



(6,0,0,0,2,0). 



(114) 



Upon decomposing it, one finds 

*„^a e R = E e A3 © V e v. 



(115) 



Two additional degrees of freedom appear in this decom- 
position when compared to the case of graphene {i.e. the 
scalar E and the pseudoscalar A3). As we explain be- 
low, they give rise to high energy bands and are naturally 
projected out for generic values of the couplings in the 
resulting effective Lagrangian. 

Quadratic terms : 

We now turn to the construction of the possible low 
energy effective terms in the Lagrangian. Having defined 
the transformation properties of "^ma under De, these 
terms are obtained from the projection of the appro- 
priate field products into De singlets. As discussed for 
the graphene case, crystalline momentum conservation is 
equivalent to demanding that the terms entering into the 
Lagrangian must not mix the Dirac points. 

The decomposition of a fermion bilinear results in 



^Ina^ ma ^ R ^ 



)R = 4Ee2Aie2A2©4A3©6V( 



)6V'. 
(116) 

Following the same steps as in the case of graphene, the 
resulting effective Lagrangian to quadratic order in the 
fields reads 

£2 = /i ^"1"^ + A ^"'"M^ ^ i^^^^v ^ V ^ V 7i)9i^ , 

(117) 



where 



M 



( 

cj* 



UJ LjJ 
UJ 





UJ 



\ 



UJ 

cj* 
cj* / 



(118) 



and the explicit expression of the three sets of 7-matrices 
is given in appendix [B] 



The coefficient /i in (117) amounts to a chemical po- 



tential for the electrons. The interesting new ingredient 
in the Kagome lattice is the M invariant tensor. In the 
abasence of M (A = 0) we obtain, at zero momentum, six 
degenerate states. The existence of a non-zero A allows 
for a splitting among this 6 states into 4 + 2 degenerate 
states. This fact leads to the interpretation of the coeffi- 
cient A as setting a gap between the six states contained 
in ^. Therefore the degeneracy of the lowest energy state 
at zero momentum can be set to 6, 2 or 4 depending on 
whether A = 0,A>OorA<0. The lowest order correc- 
tion in momentum to the ground state is non-zero for the 
four degenerate modes (which get split into 2+2), while 
to linear order the two splitted states get no correction. 

We can now see a beautiful agreement of our group 
theory approach with the phenomenology of the Kagome 
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system: a tight binding approach apphed to fermions 
hopping to nearest neighbors in the kagome lattice [2 13] 
shows that at fihing 1/3 only four degrees of freedom ap- 
pear at low energ}£j The existence of M from the group 
theory approach allows a nice match with this result if 
we set A to be positive and very large. Notice that in 
that case the assumption that the high-energy bands can 
be expanded around the minima of the low-energy bands 
is innocuous. 

At this point, a consistent low energy expansion de- 
mands to project out the massive modes. Implementing 
this projection amounts to eli mina te the irreducible rep- 
resentations E and A3 from (1151^^, which then implies 



that the low energy effective action coincides with the 
one found in the graphene case. 

Electron-phonon coupling: 

From this simple analysis we conclude that at low en- 
ergy, the electron-phonon coupling terms for the kagome 
lattice have the same expressions as in the graphene 
case, without the need of going through detailed tight- 



binding computations. More explicitely, for large A ( 115 



matches exactly (68) and hence the decomposition (185 



would be the same, implying that the interaction terms 
would be given by ( 86 ) , ( 90 ) . 



This result shows the power of our symmetry based 
approach to derive low energy effective actions. 



IV. CONCLUSION 

In the present paper we provide a method to construct 
the effective low energy Lagrangian of an electronic sys- 
tem defined on a given lattice. We use a symmetry based 
approach that incorporates the notion of memory ten- 
sors, that are nothing but the invariant tensors of the 
point group. More specifically, we have revisited in de- 
tail how symmetries constrain the low energy dynamics 
of a generic lattice system and how to retain the neces- 
sary specific information about the underlying discrete 
symmetries. 

The main idea is that in the infrared, the discrete 
translations of the lattice become continuous transla- 
tions, while the point group retains its discrete nature. 
In other words, the low energy dynamics is described by 
a field theory with fields transforming in representations 
of the discrete point group Gp. We have shown how to 
construct systematically the infrared effective Hamilto- 
nian, invariant under the full symmetry group, encoding 



the discrete origin of the system into the invariant ten- 
sors. In other words, the low energy dynamics of a sys- 
tem of electrons in a given lattice would be determined 
by the content of memory tensors of Gp and that allows 
to construct all the different invariant local terms in the 
effective Hamiltonian. 

We have applied the method developed to the case of 
the De point group, and focused into both graphene and 
kagome systems. For the graphene case we have shown 
how to derive the low energy effective field theory re- 
producing the tight-binding expressions for the electron- 
phonon system [12] jl^. In a second step, we have ap- 
plied the method to electrons in the kagome lattice close 
to filling 1/3 mis]. In this example, the method shows 
its power since the derivation of the effective action is 
quite straightforward. 

The advantages of the present approach is that it can 
be applied to study more general cases, i.e. lattices with 
other point groups, higher order in the fields and/or their 
derivatives, in a straightforward manner. 

Using the symmetry based approach presented in this 
paper one could envisage the study of situations in which 
lattice deformations contribute to next to leading orders 
in a systematic way. We believe that, for the particular 
case of two-dimensional systems, as graphene, it is possi- 
ble to overcome the small deformations expansion. These 
ideas will be explored in future investigations. This 
could be a particularly important step forward in mod- 
eling some experimental setups in which deformation of 
graphene sheets goes beyond small perturbations of the 
fiat lattice. Another possible generalization that can be 
easily envisaged within the present symmetry based set 
up is that of the inclusion of boundaries and point-like 
defects. Our procedure could also be applied to three di- 
mensional systems, such as hyper- pyrochlore lattice stud- 
ied in 3 . 
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^ The dispersion relation for the three bands found in [2| ^ show 
two bands similar to the ones found in graphene and a third one 
being flat and separated from the flrst two by a gap set by the 
hopping parameter. 

^ One can check that the two degenerate modes of M correspond 
to the E and A3 inside ^. 
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Appendix A: An alternative approach for 
determining the fermion transformation properties 

In this appendix we show that an alternative approach 
could be used to derive the transformation properties of 
the low energy degrees of freedom. 

Graphene: From the hypothesis of two inequivalent Dirac 
points, we expect four degrees of freedom at low energy. 
They correspond to: (i) the two atoms inside the funda- 
mental cell which we denote n = a, 6, and (ii) the two 
inequivalent Dirac points in Fourier space at which the 
Fermi surface degenerates, which we denote N = 1,2 
(necessarily at the vertices of the BZ). With the addi- 
tional input of a linear dispersion relation for the low 
energy degrees of freedom, we will be able to deduce the 
same transformation properties as obtained in the text. 

We start by constructing a linear representation of the 
De point symmetry group acting on these four degrees of 
freedom. To this end, we place them inside a 4-tuple ^ 



(Al) 



To fix the transformation properties of these degrees of 
freedom under Gp we orient the hexagon so that two 
of its vertices lay on the y-ax.is. When Fourier trans- 
forming, the first Brillouin zone consists of an hexagon, 
now rotated 90^. Our conventions coincide with those in 
[21 ,[23 . 

It is immediate to see that a fundamental rotation R 
interchanges both the atoms and the Dirac points. 



a ^ b, 1^2 
We therefore propose 



Rai 
Rb2 



- ab2, 
: f3ai , 



(A2) 



(A3) 



with phases a, /3 to be determined below from the point 
group constraints in (33) and the linear dispersion rela- 
tion conditions. 

Time reversal T symmetry fixes the transformations 
of the remaining degrees of freedom (a2,6i) as we now 
explain. First notice that in the absence of spin, T sim- 
ply acts as complex conjugation. Taking into account the 
dependence of the wave function on momentum one con- 
cludes that time reversal interchanges the Dirac points 
and changes the sign of the momentum with no action 



on the atoms positions (see [21] for a detailed discussion) . 
One therefore has 



Tbi=b2. 



Acting with T on (A3), one obtains 



Ra2 = a* bi 
Rbi =/3*a2. 



(A4) 



(A5) 



The phases a, /3 are not completely arbitrary, we should 
demand that the transformation D{R) defined by (A3)- 
(A5) satisfies the group constraints ([33]), that is 



{D{R)f^ = ^ 



1 



(A6) 



Let us now analyze the action of reflection on 

the fermions. For the cell orientation we have chosen, 
it is easy to see that it preserves the Dirac points but 
interchanges the atoms. 



Fx ai bi 
Fx bi= ai. 



(A7) 
(A8) 



The transformation of (a2,62) is again fixed by time re- 
versal to be 



Fx a2 = b2 
Fx ^2 = ^2 . 



(A9) 
(AlO) 



The definitions (A8)-(A10) automatically satisfy the sec- 
ond constraint in (|33D 



(All) 



(Z)(F,))2* = *. 



The remaining non-trivial constraint among the De gen- 
erators 



FxRFx — R 



(A12) 



gives no new constraint on the phases other than the one 
we have already found in (A6). 



De on ^, 



From ( |A6| ) one finds three possibilities for representing 
a(3 = 1, (A13) 



or 



af3 = e 



(A14) 



These last two possibilities can be shown to be equivalent 
upon computing the character of the representatiorj^ 

Now the requisite of a linear dispersion relation at low 
energy for the fermions comes into play. This condition 



■'^■'^ In what follows we will omit the momentum dependence of ^. 



12 The character for D is x = (4, 0, 0, 2(a/5 + 0, 0)). 
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requires an effective action containing a De singlet con- 
structed out of two fermions and one spatial derivative. 
Since the derivative transforms in the representation V, 
( [41] ) implies that the only possibility for a singlet comes 
from contracting d with another V representation. The 
second choice in (A14) guarantees that this will be pos- 



sible (in principle two possible singlets can be obtained 
when combining two fermions with a derivative see the 
textj^ Hence, the fermion dispersion relation at low 
energy dictates the choice of phases for the representa- 
tion^^ In a generic case, a linear dispersion relation for 
fermions would be forbidden by De symmetry, if no pos- 
sible choice of phases could lead to a V component in the 
fermions bilinear decomposition. 



Summarizing: the fermionic degrees of freedom (Al) mix 
under the point symmetry group as 



(A17) 



with 



D{R) 



f a 

a 



a* 



V 



(A18) 



( 
1 

VI 



/ 

i\ 

1 



(A19) 



where a = e^3 corresponds to a 60^ rotation. The repre- 



sentation (A18)-(A19) is unitary, reducible and decom- 
poses in De irreducible representations as [22 (cf. ([68|) 

^ = V © V (A20) 



Kagome: We now apply the procedure above to the 
kagome lattice. The kagome lattice has three atoms a, 6, c 
per fundamental cell (n = 3), and as for the graphene 
case we consider the simplest non-trivial case of a (degen- 
erate) band structure with two inequivalent Dirac points 



(A^ = 2). We place the six degrees of freedom inside a 
6-tuple as 



(A21) 



ci 
^1 
a2 

Under a fundamental rotation R the fermions mix simul- 
taneously as 



a ^ c ^ a, 1^2. 
Allowing for possible phases, this amount to 



Ka\ = ab2 
Rb2 =f3ci 
Rci =7a2. 



r 



Ra2 = a* bi 
Rbi =/3*C2 
Rc2 =7*^1, 



(A22) 



(A23) 



where as before the phases will be determined from the 
constraints (33). A reflection through the x-axis pre- 



serves the cones, leaves the c degrees of freedom invariant 
and permutes a ^ b^ that is 



Fx ai bi 
Fx bi ai 
Fx ci = ci . 



(A24) 
(A25) 
(A26) 



The remaining degrees of freedom transformation prop- 
erties are fixed by time reversal and amount to similar 
expressions changing 1^2. 



The first two constraints in ( 33 ) are automatically sat- 
isfied, and the constraint ( A12) results in a undetermined 
and = 1. Computing the character for the repre- 
sentation one finds that it is in fact independent of the 
product /3^. It is easily seen that this is due to the fact 
that all phases can be eliminated by appropriate defini- 
tions of the degrees of freedom. The simplest choice is 
a = /3 = 7 = 1. 

Upon computing the character one finds that the 
fermionic degrees of freedom (A21) provide a representa- 



tion for De which coincides with the one deduced from the 



transformation properties of the Wannier states (115). 



We realize the appropriate choice ( |A14| > as 

a = i3 = uj = e'^. (A15) 

A fortiori, one can check that the representation ( |A18| ), |A19| ) 
is equivalent to that proposed in [21 . Had we chosen the first 
possibility in ( |A13| the decomposition of ^ would have been 

^ G E e Al e A2 e A3 . (ai6) 

which does not allow for a standard kinetic term for the fermions 
in the sense explained above 
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Appendix B: Kagome gamma matrices 

First set: 

/ -1/4 \ 



71 



-1/4 



1/2 



1/4 



1/4 



-1/2/ 



/ V3/4 



72 



-73/4 



-V3/4 



V3/4 



0/ 



(Bl) 



(B2) 



Second set: 
/ 



71 



-e-'3 1 
e"^i e"^¥ 
1 e^^ 



/ 



72 ^ 



6 6 6 2 



6-^6 



e'3 
e-'t -1 
e-^t -10/ 

\ 



(B3) 



(B4) 



Third set: 



71 



/ V3e^^ e^t 



/ 

1 



72 



2 ^ 2^ 



" 2 2^3 





1 + 



1 ^ 5i 
"2 ^ 2^ 



V3 2 2^ 



^e^t 



"2 ^ 2^ 2 ^ 2^ 



2 2^3 ^ 

1 _ 52 -1 _ _22_ 

"2 2^ V3 



1 + 



V3 

/ 



(B5) 



(B6) 
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